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We  demonstrate  a  novel  artificial  optical  material,  the  “photonic  hyper  crystal”,  which  combines  the  most 
interesting  features  of  hyperbolic  metamaterials  and  photonic  crystals.  Similar  to  hyperbolic  metamaterials, 
photonic  hyper  crystals  exhibit  broadband  divergence  in  their  photonic  density  of  states  due  to  the  lack  of 
usual  diffraction  limit  on  the  photon  wave  vector.  On  the  other  hand,  similar  to  photonic  crystals,  hyperbolic 
dispersion  law  of  extraordinary  photons  is  modulated  by  forbidden  gaps  near  the  boundaries  of  photonic 
BrUlouin  zones.  Three  dimensional  self  assembly  of  photonic  hyper  crystals  has  been  achieved  by 
application  of  external  magnetic  field  to  a  cobalt  nanoparticle  based  ferrofluid.  Unique  spectral  properties 
of  photonic  hyper  crystals  lead  to  extreme  sensitivity  of  the  material  to  monolayer  coatings  of  cobalt 
nanoparticles,  which  should  find  numerous  applications  in  biological  and  chemical  sensing. 


Over  the  last  few  decades  a  considerable  progress  has  been  made  in  developing  artificial  optical  materials 
with  novel  and  often  counterintuitive  properties.  Revolutionary  research  by  Yablonovitch  and  John  on 
photonic  crystals'  "  was  followed  by  the  development  of  electromagnetic  metamaterial  paradigm  by 
Pendry".  Even  though  considerable  difficulties  still  exist  in  fabrication  of  three  dimensional  (3D)  photonic 
crystals  and  metamaterials,  both  fields  exhibit  considerable  experimental  progress'''".  On  the  other  hand,  on 
the  theoretical  side  these  fields  are  believed  to  be  complementary  but  mutually  exclusive.  Photonic  crystal  effects 
typically  occur  in  artificial  optical  media  which  are  periodically  structured  on  the  scale  of  free  space  light 
wavelength  X,  while  electromagnetic  metamaterials  are  required  to  be  structured  (not  necessarily  in  a  periodic 
fashion)  on  the  scale,  which  is  much  smaller  than  the  free  space  wavelength  of  light.  For  example,  in  metal 
nanowire  based  hyperbolic  metamaterials'’  schematically  shown  in  Fig.  1 A  the  inter  wire  distance  must  be  much 
smaller  than  /,.  Here  we  report  experimental  realization  of  3D  “photonic  hyper  crystals”  which  bridge  this  divide 
by  combining  the  most  interesting  properties  of  hyperbolic  metamaterials  and  photonic  crystals. 

Our  concept  of  the  photonic  hyper  crystal"  is  based  on  the  fact  that  dispersion  law  of  extraordinary  photons  in 
hyperbolic  metamaterials" 


kl  ,  kl  +  kj 


(1) 


does  not  exhibit  the  usual  diffraction  limit.  In  such  uniaxial  metamaterials  the  in  plane  8xy  and  out  of  plane  Sz 
components  of  the  dielectric  permittivity  tensor  have  opposite  signs  (e.g.  the  metal  wire  array  hyperbolic  meta 
material  shown  in  Fig.  1 A  may  have  Sz  <  0  and  >  0  "),  so  that  the  photon  wave  vector  components  fc,  are  not 
bounded  at  a  given  frequency  of  light  ro.  Existence  of  large  k  vector  modes  in  a  broad  range  of  frequencies  is 
responsible  for  such  unusual  effects  as  hyperlens  based  super  resolution  imaging®  "  and  broadband  divergence 
of  photonic  density  of  states  in  hyperbolic  metamaterials'".  On  the  other  hand,  this  also  means  that  periodic 
modulation  of  hyperbolic  metamaterial  properties  on  a  scale  L  <^X  (see  inset  in  Eig.  1  A)  would  lead  to  Bragg 
scattering  of  extraordinary  photons  and  formation  of  photonic  band  structure  no  matter  how  small  L  is®.  Thus,  so 
formed  “photonic  hyper  crystals”  would  combine  the  most  interesting  properties  of  hyperbolic  metamaterials 
and  photonic  crystals.  For  example,  similar  to  classic  photonic  crystal  effect  predicted  by  John",  strong  local 
ization  of  photons  may  occur  in  photonic  hyper  crystals.  However,  unlike  usual  photonic  crystals  where  light 
localization  occurs  on  a  scale  ~  7,  photonic  hyper  crystals  may  exhibit  light  localization  on  deep  subwavelength 
scale.  Similar  to  surface  plasmon  resonance  (SPR)'"  and  surface  enhanced  Raman  (SERS)''*  based  sensing, 
engineered  localization  of  light  on  deep  subwavelength  scale  in  photonic  hyper  crystals  should  find  numerous 
applications  in  biological  and  chemical  sensing. 
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Figure  1  |  (A)  Experimental  geometry  of  the  ferrofluid  based  hyperbolic  metamaterial.  The  array  of  self  assembled  cobalt  nanocolumns  has  typical 
separation  a~  20  nm  between  the  nanocolumns.  The  inset  shows  a  photonic  hyper  crystal  structure  formed  by  periodic  arrangement  of  cobalt  rich  and 
cobalt  sparse  regions  with  typical  periodicity  L~  2  pm,  so  that  periodic  modulation  of  hyperbolic  metamaterial  properties  on  a  scale  L  ^  X  is  achieved  in 
the  LWIR  spectral  range  where  2  ~  10  pm.  Since  photon  wave  vector  in  hyperbolic  metamaterials  is  not  diffraction  limited,  periodic  modulation  of 
hyperbolic  metamaterial  properties  on  a  scale  L  ^  X  would  lead  to  Bragg  scattering  and  formation  ofband  structure.  (B  F)  Microscopic  images  of  cobalt 
nanoparticle  based  ferrofluid  reveal  subwavelength  modulation  of  its  spatial  properties:  frames  (B)  and  (D)  show  microscopic  images  of  the  diluted 
cobalt  nanoparticle  based  ferrofluid  before  and  after  application  of  external  magnetic  field.  The  pattern  of  self  assembled  stripes  visible  in  image  D  is  due 
to  phase  separation  of  the  ferrofluid  into  cobalt  reach  and  cobalt  poor  phases.  The  stripes  are  oriented  along  the  direction  of  magnetic  field.  The  inset 
shows  Fourier  transform  image  of  frame  D.  Its  cross  section  presented  in  panel  E  shows  a  histogram  of  different  periods  present  in  the  image.  A  movie 
taken  during  this  experiment  when  the  external  magnetic  field  was  turned  on  and  off  is  available  in  the  supporting  on  line  material.  A  microscopic  image 
of  the  sample  taken  along  the  axes  of  the  nanowires  is  shown  in  frame  C.  Panel  F  demonstrates  that  the  original  undiluted  ferrofluid  exhibits  similar  phase 
separation  in  external  magnetic  field,  even  though  on  much  smaller  scale. 
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Band  structure  and  field  distribution  inside  a  photonic  hyper 
crystal  may  be  obtained  in  a  straightforward  manner.  While  both 
and  may  exhibit  periodic  spatial  dependencies,  let  us  consider 
the  relatively  simple  case  of  coordinate  independent  e^y  >  0  and 
periodic  e^{z)  <  0  with  a  period  L  <g  1.  Aside  from  the  relative 
mathematical  simplicity  of  this  model,  it  also  corresponds  to  the 
most  readily  available  low  loss  realizations  of  a  hyperbolic  metama 
terials  such  as  the  composites  formed  by  metallic  nanowires  in  a 
dielectric  membrane'^  (where  e^y  >  0  and  e^{z)  <  0),  and  planar 
layered  metal  dielectric'’’  '’’  and  doped  semiconductor  metamater 
ials'®.  Taking  into  account  the  translational  symmetry  of  the  system 
in  X  and  y  directions,  we  can  introduce  the  in  plane  wave  vector 
ky)  so  that  the  propagating  waves  can  be  expressed  as 

Eoj{r)  =  E(z)  exp  (ik^x  +  ikyy)  (2) 

Dco(r)  =  D(z)  exp  (ik^x  +  ikyy) 

Bai(r)  =  B(z)  exp  +  ikyy) 


The  uniaxial  symmetry  of  this  medium  reduces  the  ordinary  and 
extraordinary  waves  to  respectively  the  TE  (EEz)  and  TM  (B-Lz) 
polarized  modes.  Introducing  the  wavefunction  as  the  z 

component  of  the  electric  displacement  field  of  the  TM  wave 

il/(r)  =  D^(r)  =  £jz)E^(r)=  - -k^B  (3) 

CO 


for  the  wave  equation  we  obtain 


^  =  £xy 


(4) 


In  this  wave  equation  the  periodic  ratio  acts  as  a  periodic 
effective  potential.  As  usual,  solutions  of  eq.(4)  may  be  found  as 
Bloch  waves 


00 

i/'(z)=  ^mexp{iih  + 

m  0 


(5) 


where  is  defined  within  the  first  BrUlouin  zone  —K/L<kz<n/L. 
Strong  Bragg  scattering  is  observed  near  the  Brillouin  zone  bound 
aries  at  k^  ~  n/L  ^  kI).,  leading  to  the  formation  of  photonic  band 
gaps  in  both  the  wavenumber  and  the  frequency  domains.  This 
behavior  is  illustrated  in  Fig.  2,  where  we  compare  the  dispersion 
diagram  for  an  example  of  such  nanowire  based  photonic  hyper 
crystal  to  its  effective  medium  counterpart.  The  material  parameters 
of  photonic  hyper  crystal  are  based  on  the  parameters  of  stratified 
ferrofluid  described  below.  However,  we  should  emphasize  that  Fig. 
2A  is  presented  as  an  illustration  only,  since  in  real  experiment  we 
cannot  separate  magnetic  field  induced  hyperbolic  properties  of  the 
ferrofluid  from  its  photonic  hyper  crystal  properties  (or  vice  versa), 
at  the  actual  wavelengths  where  both  properties  coexist. 

Similar  to  the  usual  photonic  crystals^,  adiabatic  chirping  of  L 
leads  to  strong  field  enhancement  which,  unlike  that  in  the  conven 
tional  photonic  crystals,  occurs  on  a  deep  subwavelength  scale.  We 
should  also  mention  that  the  most  interesting  case  appears  to  be  the 
epsilon  near  zero  (ENZ)  situation  where  approaches  zero  near  a 
periodic  set  of  planes.  As  has  been  demonstrated  in  Ref  19,  electric 
field  of  the  extraordinary  wave  diverges  in  these  regions.  These  peri 
odic  field  divergences  appear  to  be  most  beneficial  for  sensing 
applications. 

In  order  to  validate  the  photonic  hyper  crystal  concept  we  have 
developed  an  experimental  technique  which  uses  three  dimensional 
self  assembly  of  cobalt  nanoparticles  in  the  presence  of  external 
magnetic  field.  Magnetic  nanoparticles  in  a  ferrofluid  are  known  to 
form  nanocolumns  aligned  along  the  magnetic  field^".  Moreover, 


depending  on  the  magnitude  of  magnetic  field,  nanoparticle  concen 
tration  and  solvent  used,  phase  separation  into  nanoparticle  rich  and 
nanoparticle  poor  phases  may  occur  in  many  ferrofluids^'.  This 
phase  separation  occurs  on  a  0.1  1  micrometer  scale.  Therefore,  it 
can  be  used  to  fabricate  a  self  assembled  photonic  hypercrystal. 

For  our  experiments  we  have  chosen  cobalt  magnetic  fluid  27 
0001  from  Strem  Chemicals  composed  of  10  nm  cobalt  nanoparti 
cles  in  kerosene  coated  with  sodium  dioctylsulfosuccinate  and  a 
monolayer  of  LP4  fatty  acid  condensation  polymer.  The  average 
volume  fraction  of  cobalt  nanoparticles  in  this  ferrofluid  is  p 
8.2%.  Cobalt  behaves  as  an  excellent  metal  in  the  long  wavelength 
infrared  range  (LWIR),  as  evident  hy  Fig.  3A:  the  real  part  of  its 
refractive  index,  n,  is  much  smaller  than  its  imaginary  part, 
Thus,  real  part  of  e,  Ree  rE  k^,  is  negative,  and  its  absolute  value 
is  much  larger  than  its  imaginary  part,  Ime  2nk.  Therefore,  it  is 
highly  suitable  for  fabrication  of  hyperbolic  metamaterials.  The 
structural  parameter  of  such  a  metamaterial  falls  into  a  few  nan 
ometers  range;  the  cobalt  nanoparticle  size  is  10  nm,  while  average 
inter  particle  distance  at  8.2%  volume  fraction  is  about  19  nm. 
Therefore,  the  metamaterial  properties  may  be  described  by  effective 
medium  parameters  on  spatial  scales  —100  nm.  On  the  other  hand 
as  demonstrated  below,  ferrofluid  begins  to  exhibit  hyperbolic  beha 
vior  in  the  range  of  free  space  wavelengths  —10000  nm  and  above 
in  the  so  called  long  wavelength  infrared  (LWIR)  frequency  range. 
Thus,  in  between  100  nm  and  10000  nm  there  exists  an  ample  range 
of  spatial  scales  which  enable  photonic  hyper  crystal  behavior 
described  above.  For  example,  if  the  effective  medium  parameters 
of  ferrofluid  are  modulated  on  the  scale  of  —2000  nm  (as  evident 
from  Fig.  ID)  the  large  k  vector  modes  which  exist  in  the  metama 
terial  in  the  LWIR  range  will  experience  Bragg  scattering  due  to  this 
modulation. 

Electromagnetic  properties  of  these  metamaterials  may  be  under 
stood  based  on  the  Maxwell  Garnett  approximation  via  the  dielectric 
permittivities  £„  and  of  cobalt  and  kerosene,  respectively.  Volume 
fraction  of  cobalt  nanoparticles  aligned  into  nanocolumns  by 
external  magnetic  field,  a(B),  depends  on  the  field  magnitude.  At 
very  large  magnetic  fields  all  nanoparticles  are  aligned  into  nanocol 
umns,  so  that  0!(oo)  p  8.2%.  We  will  assume  that  at  smaller  fields 
the  difference  0!(oo)  a(B)  describes  cobalt  nanoparticles,  which  are 
not  aligned  and  distributed  homogeneously  inside  the  ferrofluid. 
Dielectric  polarizability  of  these  nanoparticles  may  be  included  into 
ej  leading  to  slight  increase  in  its  value.  Using  this  model,  the  diag 
onal  components  of  the  ferrofluid  permittivity  may  be  calculated 
using  Maxwell  Garnett  approximation  as  follows^^: 

£^  =  a(B)£^  +  (l-a(B))£d  (6) 

^  2v.(B)£^£d  +  (1  -  a(B))£d{£d  +  £,„) 

(l-a{B)){£d  +  £„,)  +  2a{B)£d 


Calculated  wavelength  dependences  of  £z  and  £xy  at  a(oo)  p 
8.2%  are  plotted  in  Eigs.  3CD.  While  8xy  stays  positive  and  almost 
constant,  £^  changes  sign  to  negative  around  a  3  pm.  If  the  volume 
fraction  of  cobalt  nanoparticles  varies,  this  change  of  sign  occurs  at 
some  critical  value  0!h: 

a.>a.H= -  (8) 

e-d-e-m 

The  value  of  ocu  as  a  function  of  wavelength  is  plotted  in  Fig.  3B.  This 
plot  indicates  that  the  original  ferrofluid  diluted  with  kerosene  at  a 
1:10  ratio  remains  a  hyperbolic  medium  above  a  5  pm.  More 
interestingly,  such  a  diluted  ferrofluid  develops  very  pronounced 
phase  separation  into  periodically  aligned  cobalt  rich  and  cobalt  poor 
phases  (with  periodicity  L  —  2  pm)  if  subjected  to  external  magnetic 
field.  Optical  microscope  images  of  the  diluted  ferrofluid  before  and 
after  application  of  external  magnetic  field  are  shown  in  Figs.  IB  and 
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Figure  2  |  Comparison  of  the  effective  medium  dispersion  of  a  nanowire  based  hyperbolic  metamaterial  (A)  to  the  exact  solution  for  photonic  crystal 
(B).  The  hyper  crystal  unit  cell  in  (B)  is  assumed  to  be  1000  nm,  with  one  half  of  the  cell  is  filled  with  the  same  ferrofluid  while  another  half  is  pure 
kerosene.  The  ferrofluid  aligned  by  external  magnetic  field  applied  in  the  z  direction,  is  characterized  by  the  material  parameters  described  in  Fig.  3.  The 
stratified  ferrofluid  is  assumed  to  form  layers  with  the  normal  along  the  x  direction.  The  effective  medium  dispersion  in  panel  (A)  uses  the  dielectric 
permittivity  tensor  obtained  by  the  homegenization  of  the  electromagnetic  response  of  the  hyper  crystal  unit  cell.  This  is  an  artificial  result  that  would  be 
expected  if  we  could  experimentally  isolate  magnetic  field  induced  hyperbolic  behaviour  of  the  ferrofluid  from  its  photonic  hyper  crystal  behaviour.  In 
both  panels,  the  wavevector  components  are  given  in  units  of  the  free  space  wavenumber  ko. 


ID.  The  periodic  pattern  of  self  assembled  stripes  visible  in  image  D 
appears  due  to  phase  separation.  The  stripes  are  oriented  along  the 
direction  of  magnetic  field.  A  movie  taken  during  this  experiment 
when  magnetic  field  was  turned  on  and  off  is  available  in  the  sup 
porting  on  line  material.  It  clearly  demonstrates  the  tunability  of  the 
medium  by  external  magnetic  field.  The  stripe  periodicity  L  ~  2  pm 
appears  to  be  much  smaller  than  the  free  space  wavelength  in  the 
hyperbolic  frequency  range.  Therefore,  created  self  assembled  opti 
cal  medium  appears  to  be  a  photonic  hyper  crystal.  We  should  also 
note  that  the  original  undiluted  ferrofluid  exhibits  similar  phase 
separation  in  external  magnetic  field,  even  though  on  much  smaller 
—0.3  pm  spatial  scale  (see  Fig.  IF). 

Polarization  dependencies  of  ferrofluid  transmission  as  a  function 
of  magnetic  field  and  nanoparticle  concentration  measured  in  a 
broad  0.5  pm  16  pm  wavelength  range  conclusively  prove  hyper 
bolic  crystal  character  of  ferrofluid  anisotropy  in  the  long  wavelength 
IR  range  at  large  enough  magnetic  field.  Fig.  3E  shows  polarization 
dependent  transmission  spectra  of 200  pm  thick  undiluted  ferrofluid 
sample  obtained  using  FTIR  spectrometer.  These  data  are  consistent 


with  hyperbolic  character  of  e  tensor  of  the  ferrofluid  in  B  1000  G. 
Ferrofluid  transmission  is  large  for  polarization  direction  perpendic 
ular  to  magnetic  field  (perpendicular  to  cobalt  nanoparticle  chains) 
suggesting  dielectric  character  of  £  in  this  direction.  On  the  other 
hand,  ferrofluid  transmission  falls  to  near  zero  for  polarization  dir 
ection  along  the  chains,  suggesting  metallic  character  of  £  in  this 
direction.  However,  these  measurements  are  clearly  affected  by 
numerous  ferrofluid  absorption  lines,  whose  behaviour  wHl  be 
addressed  later.  In  order  to  be  conclusive,  LWIR  data  must  be  sup 
ported  by  polarization  measurements  in  the  visible  range  presented 
in  Fig.  4.  Main  features  of  these  dependencies  may  indeed  be  under 
stood  based  on  the  Maxwell  Garnett  approximation  (eqs.(6,7)).  At 
small  magnetic  fields  ot(B)  <C  p  and  we  obtain 

£z  =  a£m  +  (l  — a)£d— fid  +  Me'm  and  Exy~£d,  (9) 

where  £"„  is  the  imaginary  part  of  £m-  This  absorption  related  aniso 
tropy  is  primarily  responsible  for  the  polarization  dependence  of 
ferrofluid  transmission.  Results  presented  in  Figs.  4AB  were  obtained 
at  2  1.55  pm  by  measuring  transmission  of  ferrofluid  placed  inside 
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Figure  3  |  (A)  Optical  properties  of  cobalt  as  tabulated  in  Ref.  22:  real  (n)  and  imaginary  (k)  parts  of  the  cobalt  refractive  index  are  plotted  in  the  long 
wavelength  IR  range.  (B)  Critical  volume  fraction  of  cobalt  nanoparticles  corresponding  to  ferrofluid  transition  to  hyperbolic  metamaterial  phase 
shown  as  a  function  of  free  space  light  wavelength.  (C,D)  Wavelength  dependencies  of  Sj.  and  at  o:(=o)  =  8.2%.  While  Ej^y  stays  positive  and  almost 
constant,  E^  changes  sign  around  X  =  3  pm.  (E)  Polarization  dependent  transmission  spectra  of  200  pm  thick  ferrofluid  sample  measured  using  FTIR 
spectrometer  are  consistent  with  hyperbolic  character  of  E  tensor. 


ad  10  pm  wide  quartz  cuvette  (note  that  fiquauz  ~  £d  in  the  visible 
and  near  IR  ranges).  If  ferrofluid  inside  the  cuvette  may  be  consid 
ered  as  a  homogeneous  medium,  its  transmission  as  a  function  of 
light  polarization  (|)  should  behave  as 

TpF  =  Tmin  sin^  +  Tmax  cos^  ^  1  sin^  (p^  (10) 

where  Tmin  is  transmission  for  polarization  direction  parallel  to  mag 
netic  field  and  T^ax  is  transmission  for  polarization  direction  per 


pendicular  to  magnetic  field.  Since  electric  field  equals  E  =  Eq  exp 
/  2n^s  \ 

(  i  — - —  d  I  after  passing  through  a  ferrofluid  layer  of  thickness  d, 

the  contrast  between  transmissions  can  be  estimated  using  expres 
sions  for  and  from  eq.(9): 


In 


2ndcc  ( 


(11) 
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Figure  4  |  (A,B)  Experimentally  measured  transmission  of  the  cobalt  based  ferrofluid  at  A,  =  1.55  pm  as  a  function  of  cobalt  concentration,  external 
magnetic  field  and  polarization  angle.  As  shown  in  panel  (B),  in  the  absence  of  phase  separation  polarization  curves  exhibit  the  Malus  law  sitf  (j) 
dependencies.  On  the  other  hand,  in  the  phase  separated  state  polarization  curves  exhibit  universal  “polarization  notch”  behaviour  indicating  long  range 
order  of  periodically  aligned  cobalt  nanocolumns.  (C)  Experimentally  measured  transmission  of  the  ferrofluid  at  a('X))  =  8.2%,  X  =  488  nm  and 
B  =  1630  G.  (D)  Polarization  contrast  measurements  (data  points)  in  the  0.5  10.6  pm  range  provide  experimental  validation  of  Maxwell  Garnett 
approximation  (red  line)  at  low  nanoparticle  concentrations.  (E)  Visible  and  near  IR  light  propagation  through  a  stratified  ferrofluid.  Due  to  much  higher 
absorption  in  the  cobalt  rich  phase,  ferrofluid  may  be  represented  as  multiple  Fabry  Perot  like  resonators  separated  by  strongly  absorptive  regions  of 
cobalt  rich  phase.  (F)  Calculated  transmission  of  such  a  model  stratified  ferrofluid  at  =  1.55  pm. 


Equation  (11)  may  be  used  to  verify  Maxwell  Garnett  approximation 
and  measure  a(B). 

At  A  1.5  pm  and  large  fields,  polarization  curves  depart  from  the 
Malus  law,  showing  much  sharper  minima  (Fig.  4A),  which  will  be 
discussed  later.  However,  in  the  visible  range,  or  in  smaller  fields, 
polarization  curves  indeed  obey  the  Malus  law  exhibiting  sin^i^ 
dependence  (see  Fig.  4C).  Moreover,  low  field  measurements  m 
Fig.  4D  performed  on  the  same  sample  at  four  different  wavelengths 
across  a  very  broad  0.5  10.6  pm  wavelength  range  foUow  the 
Maxwell  Garnett  result  given  by  Eq.  1 1  (note  that  a.  0.0004  <C  p 
in  this  case,  since  only  a  small  fraction  of  particles  are  aligned  into 
nanocolumns  in  a  week  magnetic  field).  These  measurements  give 
substantial  support  to  our  model. 

Let  us  demonstrate  that  these  measurements  conclusively  point 
towards  hyperbolic  character  of  ferrofluid  anisotropy  in  the  long 
wavelength  IR  range.  As  an  example,  let  us  calculate  ai(B))  for  the 


polarization  curve  in  Fig.  4C  measured  at  B  1630  G.  While  polar 
ization  contrast  is  quite  large  for  this  curve:  Tmax/Tmm  ~  10,  it  fits 
Malus  law  sin^^  dependence  rather  well.  Based  on  Eq.  1 1,  calculated 
oi(B))  equals  0.002  in  this  experiment  performed  at  A  488  nm. 
Therefore,  at  A  14  pm  where  Re(£in)  2100  the  real  parts  of 
the  dielectric  permittivity  tensor  components  are 

+  -  2.2,  while  Re(£j:;,)«£d~ +  2.1  (12) 

meaning  that  the  magnetized  ferrofluid  does  indeed  exhibit  hyper 
bolic  metamaterial  behaviour  on  spatial  scales  above  ~100  nm. 

On  the  other  hand,  polarization  dependent  transmission  of  phase 
separated  ferrofluid  samples  shown  in  Fig.  ID  supports  long  range 
order  of  periodically  aligned  (with  periodicity  L  ~  2  pm)  cobalt  rich 
and  cobalt  sparse  regions  in  a  magnetized  ferrofluid.  U sing  Eqn.  (11), 
experimentally  observed  contrast  between  transmission  of  the  cobalt 
rich  and  cobalt  sparse  regions  in  Figs.  ID  and  IF  may  be  used  to 
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evaluate  volume  fractions  of  cobalt  nanoparticles  and  in  these 
phases.  In  both  cases  the  ratio  a/iXj  falls  into  the  1.6  1.8  range. 
Moreover,  it  is  clear  from  Fig.  4A  that  polarization  curves  measured 
at  large  magnetic  fields  deviate  strongly  from  the  Malus  law  sin^c|) 
dependence.  As  shown  in  Fig.  4B,  polarization  dependencies  of  the 
phase  separated  ferrofluid  measured  at  different  nanoparticle  con 
centrations  may  be  scaled  to  a  universal  “polarization  notch”  curve 
described  by  two  parameters: 


1  —  a  sin^  (j) 
l  +  b  cos^  (j) 


(13) 


As  illustrated  in  Fig.  4E,  such  a  “polarization  notch”  dependence 
arise  naturally  due  to  multiple  scattering  of  light  by  periodically 
spaced  nanoparticle  rich  and  nanoparticle  poor  phases  seen  in  Fig. 
ID.  Due  to  much  higher  absorption  in  the  cobalt  rich  phase  (a^  ~ 
1.7o!s),  ferrofluid  may  be  represented  as  multiple  Fabry  Perot  like 
resonators  separated  by  strongly  absorptive  regions  of  cobalt  rich 
phase.  Let  us  consider  light  propagation  through  one  of  such 
Fabry  Perot  resonators  formed  by  a  single  cobalt  sparse  region. 
Transmission  amplitude  t  through  such  an  elementary  resonator  is 
given  by  a  coherent  sum  of  aU  the  multiple  scattering  events  in  which 
the  wave  experiences  multiple  polarization  dependent  reflections 
from  the  phase  boundaries: 

00  CO  rp 

m  0  m  0  i  K.e 


where  5  is  the  phase  accumulated  by  the  wave  upon  double  pass 
between  the  phase  boundaries,  and  T  and  R  are  the  transmission 
and  reflection  coefficients  for  a  single  inter  phase  boundary  (note 
that  T  +  R  1  when  losses  at  boundaries  are  neglected).  Resulting 
transmission  through  a  single  resonator  formed  hy  the  nanoparticle 
rich  and  nanoparticle  poor  phases  is  thus 


TFP  =  tt* 


1 +R2  —  2Rcos  S 


(15) 


Since  R  exhibits  strong  polarization  dependence  due  to  alignment  of 
cobalt  nanoparticles  by  external  magnetic  field,  the  universal  “polar 
ization  notch”  curve  given  by  eq.(13)  finds  natural  explanation. 
Within  the  scope  of  model  presented  in  Fig.  4E  the  overall  transmis 
sion  of  ferrofluid  is 


nr  T"  rnN 

J  total  —  ^  FF  ^  ppf 


(16) 


where  Tpp  is  given  by  eq.(lO)  and  N  is  the  number  of  effective  Fabry 
Perot  resonators  from  Fig.  4E.  Since  the  cuvette  size  in  experiments 
presented  in  Fig.  4  is  d  10  pm,  the  number  of  resonators  N  is  small. 
Calculated  transmission  of  such  a  model  stratified  ferrofluid  is  pre 
sented  in  Fig.  4F.  For  simplicity  we  have  assumed  that  interfaces 
between  cobalt  rich  and  cobalt  sparse  regions  behave  as  ideal  polar 
izers,  so  that  R  ~  sin^^.  However,  even  such  an  oversimplified  model 
reproduces  polarization  notch  behavior  of  stratified  ferrofluid.  Note 
that  in  order  for  this  model  to  be  valid,  polarization  state  of  light  must 
be  conserved  during  propagation,  indicating  a  well  ordered  ferro 
fluid  state.  In  principle,  such  mechanism  as  Anderson  localization 
may  also  lead  to  polarization  dependent  changes  in  transmission 
(see  for  example^"')  which  study  effects  of  polarization  on  Ander 
son  localization  in  one  dimensional  random  stacks.  However,  the 
described  effect  consists  in  strong  polarization  dependencies  near 
the  Brewster  angle  (see  for  example  Fig.  11  from^^).  Therefore,  they 
are  not  applicable  to  our  experiments,  which  have  been  conducted  at 
normal  Incidence.  Moreover,  long  range  ordering  of  our  samples  is 
clearly  demonstrated  by  Fig.  1.  Relatively  small  amount  of  disorder 
visible  in  these  images  does  not  justify  application  of  Anderson  local 
ization  theory.  Thus,  our  microscopy  results  in  Fig.  1C  E,  and 
detailed  examination  of  polarization  dependent  transmission  of 
the  ferrofluid  presented  in  Fig.  4  strongly  indicate  both  hyperbolic 


character  of  the  medium  in  the  long  wavelength  IR  range,  and  its 
well  ordered  periodicity.  These  data  conclusively  prove  that  the  mag 
netized  ferrofluid  may  indeed  be  considered  as  a  self  assembled 
photonic  hyper  crystal. 

As  summarized  in  Fig.  5,  fabricated  photonic  hyper  crystals 
exhibit  all  the  typical  features  associated  with  the  hyperbolic  meta 
materials.  For  example.  Figs.  5A  and  5B  demonstrate  modulation  of 
the  metamaterial  absorption  lines  due  to  hyperbolic  order  induced  by 
the  external  magnetic  field.  Radiation  lifetime  engineering  due  to 
broadband  divergence  of  the  photonic  density  of  states  in  hyperbolic 
metamaterials*^  has  been  predicted  in^^  and  experimentally  con 
firmed  Absorption  lines  of  kerosene  offer  natural  target 

for  testing  this  effect  in  fabricated  photonic  hyper  crystals.  Absorp 
tion  spectra  measured  using  FTIR  spectrometer  with  and  without 
external  magnetic  field  are  consistent  with  the  decrease  of  the  radi 
ation  lifetime  of  kerosene  molecules  in  the  hyperbolic  state.  Shorter 
lifetime  leads  to  decrease  in  absorption  line  amplitude,  which  is 
detected  in  Fig.  5B.  Reduced  reflection  has  been  reported  recently 
as  another  experimental  signature  of  hyperbolic  metamaterials^^. 
Due  to  broadband  divergence  of  photonic  density  of  states,  rough 
ened  surface  of  a  hyperbolic  metamaterial  scatters  light  preferentially 
inside  the  medium,  resulting  in  reduced  reflectance.  In  case  of  photo 
nic  hyper  crystals,  periodic  modulation  of  the  hyper  crystal  surface 
plays  the  same  role,  facilitating  coupling  of  external  light  into  the 
large  k  vector  modes  of  photonic  hyper  crystal.  Figs.  5C  and  5D 
demonstrate  such  magnetic  field  induced  reduced  reflectivity  of 
the  ferrofluid.  Moreover,  Fig.  5D  shows  evidence  of  weak  periodic 
modulation  of  reflectivity  at  a  >  9  pm  with  period  close  to  1  micro 
meter:  reflectivity  maxima  in  Fig.  5D  are  plotted  as  a  function  of  free 
space  light  wavelength  in  Fig.  5E.  These  reflectivity  maxima  are 
consistent  with  Bragg  scattering  of  incoming  LWIR  light  due  to 
periodic  modulation  of  ferrofluid  concentration  observed  in  Fig.  1. 
Such  a  Bragg  like  scattering  is  difficult  to  explain  without  consider 
ing  large  k  vector  modes  in  hyperbolic  metamaterials.  Coinciden 
tally,  our  numerical  calculations  (see  eq.(12))  indeed  indicate 
hyperbolic  behavior  of  ferrofluid  in  the  LWIR  range. 

Finally,  we  would  like  to  illustrate  photonic  hyper  crystal  potential 
in  chemical  and  biological  sensing,  which  is  made  possible  by  spa 
tially  selective  field  enhancement  effects  described  above.  This  poten 
tial  is  revealed  by  detailed  measurements  of  magnetic  field  induced 
transmission  of  photonic  hyper  crystals  in  the  broad  IR  spectral 
range  presented  in  Figs.  6A  and  6B.  FTIR  spectral  measurements 
are  broadly  accepted  as  a  powerful  “chemical  fingerprinting”  tool 
in  chemical  and  biosensing.  Therefore,  broadly  available  magnetic 
field  tunable  photonic  hyper  crystals  operating  in  the  IR  range  open 
up  new  valuable  opportunities  in  chemical  analysis.  Our  experi 
mental  data  presented  in  Fig.  6  clearly  illustrate  this  point.  The 
FTIR  transmission  spectrum  of  the  diluted  (p  0.8%)  ferrofluid 
in  Fig.  6A  exhibits  clear  set  of  kerosene  absorption  lines,  which  is 
consistent  with  other  published  data  (see  for  example  ref  28).  On  the 
other  hand,  magnetic  field  induced  transmission  spectrum  of  the  p 
8.2%  ferrofluid  shown  in  Fig.  6B  contains  a  very  pronounced  absorp 
tion  line  at  ~  12  pm  (—840  cm  '),  which  cannot  be  attributed  to 
kerosene.  Quite  naturally  this  absorption  line  may  be  attributed  to 
fatty  acids,  since  cobalt  nanoparticles  are  coated  with  a  monolayer  of 
surfactant  composed  of  various  fatty  acids,  such  as  lactic,  oleic  etc. 
acids  as  shown  in  Fig.  6D.  Detailed  comparison  of  Fig.  6B  with  the 
10  14  pm  portion  of  lactic  acid  FTIR  absorption  spectrum  shown  in 
Fig.  6C  indeed  indicates  a  close  match.  The  fatty  acid  line  appears  to 
be  about  as  strong  as  A,  —  14  pm  (—695  cm  ')  line  of  kerosene,  even 
though  the  oscillator  strength  of  these  molecular  lines  is  about  the 
same,  while  the  amount  of  kerosene  in  the  sample  is  —  two  orders  of 
magnitude  larger  (a  monolayer  coating  of  fatty  acids  on  a  10  nm 
cobalt  nanoparticle  occupies  no  more  than  1%  of  ferrofluid  volume). 
This  paradoxical  situation  clearly  indicates  local  field  enhancement 
effects.  Another  strong  evidence  of  field  enhancement  is  provided  by 
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Figure  5  |  (A)  FTIR  transmission  spectrum  of  the  diluted  (a(“)  =  0.8%)  ferrofluid  in  zero  magnetic  field  reveals  strong  kerosene  absorption  lines. 
(B)  Measured  differential  FTIR  transmission  spectrum  T(H)  T(0)  indicate  a  decrease  in  absorption  line  amplitude,  which  is  consistent  with  the  decrease 
of  the  kerosene  radiation  lifetime  in  the  hyperbolic  metamaterial  state.  (C)  FTIR  reflection  spectrum  of  the  ferrofluid  in  zero  magnetic  field.  (D)  Measured 
differential  FTIR  reflection  spectrum  R(H)  R(0)  indicates  reduced  reflectivity  in  the  hyperbolic  state.  (E)  Reflectivity  maxima  in  Fig.  5D  plotted  as  a 
function  of  free  space  light  wavelength  show  evidence  of  Bragg  like  scattering  of  large  k  vector  hyperbolic  metamaterial  modes. 


measurements  of  extinction  coefficient  of  the  ferrofluid  presented 
in  Fig.  6E.  Ferrofluid  subjected  to  magnetic  field  exhibits  pro 
nounced  resonances  around  the  fatty  acid  absorption  line  at  X 
~  12  pm  and  the  kerosene  absorption  line  at  X  ~  14  pm. 
These  resonances  provide  clear  evidence  of  field  enhancement 


by  cobalt  nanoparticle  chains.  We  expect  that  further  optimization 
of  photonic  hyper  crystals  geometry  will  lead  to  much  stronger 
sensitivity  of  their  optical  properties  to  chemical  and  biological 
inclusions,  indicating  a  very  strong  potential  of  photonic  hyper 
crystals  in  biological  and  chemical  sensing. 
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Figure  6 1  (A)  FTIR  transmission  spectrum  of  diluted  («(<»)  =  0.8%) 
ferrofluid  exhibits  clear  set  of  kerosene  absorption  lines.  (B)  Transmission 
spectra  of  the  a(<»)  =  8.2%  ferrofluid  measured  with  and  without 
application  of  external  magnetic  field.  Magnetic  field  induced 
transmission  spectrum  contains  a  very  pronounced  absorption  line  at 
1  ~  12  pm  (~840  cm'’),  which  can  be  attributed  to  lactic  acid.  Kerosene 
absorption  lines  are  marked  with  yellow  boxes,  while  the  fatty  acid  line  at 
840  cm*’  is  marked  with  a  green  box.  (C)  The  10  14  pm  portion  of  lactic 
acid  FTIR  absorption  spectrum.  (D)  Schematic  view  of  cobalt  nanoparticle 
coated  with  a  monolayer  of  fatty  acids,  such  as  lactic  and  oleic  acid. 

(E)  Extinction  coefficient  (the  ratio  T(0)/r(90)  of  ferrofluid  transmissions 
at  0  and  90  degrees)  of  ferrofluid  subjected  to  external  magnetic  field 
exhibits  pronounced  resonances  around  the  fatty  acid  absorption  line  at  k 
~  12  pm  and  the  kerosene  absorption  line  at  14  tun.  These  resonances 

provide  clear  evidence  of  field  enhancement  by  cobalt  nanoparticle  chains. 
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